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Abstract 

Background: Mate preference behavior is an essential first step in sexual selection and is a critical determinant in 
evolutionary biology. Previously an environmental compound (the fungicide vinclozolin) was found to promote the 
epigenetic transgenerational inheritance of an altered sperm epigenome and modified mate preference 
characteristics for three generations after exposure of a gestating female. 

Results: The current study investigated gene networks involved in various regions of the brain that correlated with 
the altered mate preference behavior in the male and female. Statistically significant correlations of gene clusters 
and modules were identified to associate with specific mate preference behaviors. This novel systems biology 
approach identified gene networks (bionetworks) involved in sex-specific mate preference behavior. Observations 
demonstrate the ability of environmental factors to promote the epigenetic transgenerational inheritance of this 
altered evolutionary biology determinant. 

Conclusions: Combined observations elucidate the potential molecular control of mate preference behavior and 
suggests environmental epigenetics can have a role in evolutionary biology. 
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Background 

The current molecular paradigm for neo-Darwinian evo- 
lutionary biology is that random DNA sequence muta- 
tions, gene flow, and genetic drift promote phenotype 
variation that allows an adaptation event to facilitate nat- 
ural selection [1]. Although environment has an important 
role in natural selection, environmental factors generally 
do not have the capacity to alter DNA sequence or muta- 
tion rates. A small group of compounds can act as muta- 
gens, but the vast majority of nutritionally derived and 
environmental toxicants do not alter DNA sequence [2]. 
The current genetic paradigm does not completely explain 
many observations such as rapid evolutionary events, en- 
vironmental impacts on evolution, and the low frequency 
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of the occurrence of useful mutations [3,4]. The realization 
that epigenetics provides an additional molecular mechan- 
ism for the environment to influence genome activity and 
biology has suggested a potential role for environmental 
epigenetics in evolutionary biology [5-11]. Charles Darwin 
recognized sexual selection as one of two determinants in 
evolutionary biology, the other being natural selection [12]. 
The physical attributes and courtship rituals involved in 
mate preference are essential for reproductive fitness and 
propagation of a species. The current study examines how 
environmental factors can promote an epigenetic event to 
promote an alteration in mate preference behavior. 

Previously we demonstrated that exposure of a gestating 
female rat to an environmental compound during fetal go- 
nadal sex determination promoted epigenetic reprogram- 
ming of the male germline [13-15]. These reprogrammed 
differential DNA methylation regions (DMR) in the sperm 
epigenome have recently been shown to be induced by a 
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variety of different environmental toxicants with exposure 
specific DMR [16,17]. The initial environmental com- 
pound used was the commonly used fungicide vinclozolin 
which is an anti- androgenic endocrine disruptor [18]. The 
primordial germ cells during migration down the genital 
ridge undergo a DNA methylation erasure that then upon 
gonadal sex determination the DNA re-methylation is 
initiated in a sex-specific manner [19]. Environmental 
exposures during this developmental stage modifies the 
epigenetic programming of the male germline that be- 
comes re-programmed (imprinted-like) and promotes a 
transgenerational phenotypic variation and adult onset 
disease state in subsequent generations [13-16]. The epi- 
genetic transgenerational inheritance of adult onset disease 
(i.e. after one year of age) in males includes infertility, pros- 
tate disease, kidney disease, immune abnormalities and 
spermatogenic defects [20,21], and in females includes 
mammary tumor development, kidney disease, reproduct- 
ive tissue abnormalities and pregnancy abnormalities [22]. 
This germline mediated epigenetic transgenerational in- 
heritance of adult onset disease is mediated in part through 
alterations in the sperm epigenome [13,14]. Since the germ- 
line establishes the base line epigenome of the organism, 
all tissues in both the female and male progeny including 
the brain appear to have altered tissue specific transge- 
nerational epigenomes, transcriptomes and phenotypes 
[15,21,23-25]. 

Investigation of the epigenetic transgenerational inher- 
itance of altered brain genome activity and behaviors 
previously demonstrated anxiety-like behavior increased 
in females and decreased in males, which correlate to 
alterations in specific brain region transcriptomes [23]. 
Altered stress responses are also detected in the trans- 
generational exposure lineage animals [26]. Interest- 
ingly, previous analysis of F3 generation control and 
vinclozolin lineage female and male rats (i.e. prior to the 
onset of disease) demonstrated an alteration in mate pref- 
erence behavior [27] . The female rats, independent of con- 
trol or vinclozolin lineage, prefer control lineage males if 
given a choice. This behavioral decision raises the possi- 
bility of an epigenetic contribution to mate preference 
and sexual selection. The current study was designed to 
directly correlate the altered mate preference behavior 
with gene networks in specific brain regions in both the 
females and males. Observations elucidate the potential 
molecular control of mate preference behavior and demon- 
strates environmental factors have the capacity to promote 
the epigenetic transgenerational inheritance of altered mate 
preference. 

Systems biology analysis has allowed biological phenom- 
ena such as mate preference to be considered from the 
molecular to physiological level. The gene bionetwork 
[28] analysis previously developed to investigate the mo- 
lecular basis of disease [29] was used in the current study. 



This approach has been used to identify gene networks as- 
sociated with disease, such as obesity and diabetes [30]. 
Recently, we have used this bionetwork analysis to study a 
normal developmental process of primordial follicle devel- 
opment in the ovary [31,32]. The gene networks identified 
were found to contain growth factors that are known to 
regulate the developmental process [31,32]. These bionet- 
work analyses use a large number of microarray transcrip- 
tome analyses under different perturbations to identify 
gene clusters and modules that are coordinately regulated 
[33,34]. The gene networks observed identify the genes 
with the highest level of integration and connection (i.e. 
connectivity) that associate with the phenotype [29,33-35]. 
This genomic approach was used in the current study to 
identify the gene bionetworks in various brain regions as- 
sociated with mate preference. 

Observations demonstrate an environmental compound 
exposure can induce an epigenetic reprogramming of the 
germline that promotes epigenetic transgenerational in- 
heritance of altered mate preference behavior. Although 
no direct epigenetic modifications in the brain were exam- 
ined, the environmentally induced epigenetic transgenera- 
tional model used indicates epigenetics can be involved in 
the induction of the altered behavioral phenotypes. Sex- 
specific effects were observed in both the male and female 
brain transcriptome and behavior correlations. The gene 
networks in specific brain regions that statistically cor- 
relate with various mate preference behaviors provides 
insight into this environmentally modified transgenera- 
tional behavior. This systems biology approach has 
elucidated novel mechanisms to be considered in mate 
preference biology. 

Results 

The experimental design involved the development of 
transgenerational control and vinclozolin lineage animals 
for a mate preference behavioral analysis [27]. Subsequently, 
a transcriptome analysis was performed on 6 different brain 
regions from adult male and female F3 generation Sprague 
Dawley rats. These brain regions have previously been 
shown to be associated with mate preference behavior 
[27,36]. The transcriptome alterations were statistically 
correlated with changes in mate preference behaviors. 
As previously described [13,37], FO generation gestating 
females were transiently exposed daily to vehicle control 
DMSO or vinclozolin from embryonic day 8-14 (E8-14) 
during fetal gonadal sex determination. The Fl gener- 
ation offspring were bred at 90 days of age to generate 
F2 generation control and vinclozolin lineage progeny 
and then F2 generation animals were bred to generate 
the F3 generation control and vinclozolin lineage animals 
[13]. No sibling or cousin breeding was used to avoid any 
inbreeding artifacts. The F3 generation control and vinclo- 
zolin male and female rats were analyzed at 3-4 months 
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of age for mate preference behaviors, as previously de- 
scribed [27]. This is an age when no major adult onset 
disease has been detected or is anticipated [20]. Later at 
11 month of age, animals were sacrificed and specific 
brain regions isolated and RNA collected for microarray 
transcriptome analysis. The differentially regulated gene 
sets ("Signature lists") for each brain region were identi- 
fied. Subsequently a bioinformatics bionetwork analysis 
[23,31] was used to correlate gene modules and networks 
with mate preference behaviors observed (Figure 1). 

Females of both control and vinclozolin lineages pre- 
ferred control lineage males over vinclozolin lineage males 
[27]. The mate preference behaviors are described in detail 
in the Methods [27] (Additional file 1: Movie SI). The 
specific behavioral components associated with mate pref- 
erence analysis include the following: "Wire Mesh" in- 
volved the experimental animal investigating the stimulus 
animals directly through the Wire Mesh; "Facial Investiga- 
tion" entailed the actual nose-to-nose contact; "Plexiglas" 
refers to the experimental animal investigating the area 
immediately bordering the Wire Mesh that separated the 
experimental animal from the stimulus animal; "Walking" 
refers to general investigation of the central testing cham- 
ber as measured by undirected walking and sniffing; and 
"Still" in which the animal was stationary with minimal 
head movement. Additional file 2: Table S1A presents the 
values for each behavioral component associated with in- 
dividual animals. After completion of the behavioral ana- 
lysis the animals were sacrificed at 11 months of age and 6 
different brain regions dissected including the amygdala 
(Amy), hippocampus (Hipp), olfactory bulb (OlfB), cingu- 
late cortex (CngCtx), entorhinal cortex (EnCtx), and pre- 
optic area-anterior hypothalamus (POAH). The procedure 
to isolate the brain regions is described in the Materials 



and Methods. The isolated tissue was immediately placed 
in Trizol reagent, frozen and stored. RNA was prepared 
for microarray transcriptome analysis from each animals 
brain regions independently. 

For the microarray analysis each F3 generation control 
and vinclozolin lineage male and female animal had six 
different brain regions analyzed which totaled 134 different 
microarrays. The microarray data were pre-processed and 
demonstrated two abnormal arrays that were omitted for 
further analysis (Additional file 3: Figure SIB). Batch effect 
corrections were made for RNA preparation date and array 
scan date with no major batch effects detected. The array 
data were then processed as previously described [31] to 
identify the differentially expressed gene sets for each brain 
region (Table 1). The differentially expressed genes in the 
Signature lists required a greater than 1.2 fold change in 
expression and all changes in expression were statistically 
significant with p < 0.05, as described in the Methods. Since 
a 20% alteration in gene expression for many genes, such 
as transcription factors, can have dramatic cellular and bio- 
logical responses [26,32], a more stringent cut off (e.g. 2x) 
was not used in the current study. In the current study the 
primary focus was on the coexpression patterns of the dif- 
ferentially expressed genes through the co-expression net- 
work analysis. 

The number of control lineage versus vinclozolin lineage 
differentially expressed genes in the Signature lists ranged 
from 43 to 803 with both up-regulated and down-regulated 
genes (Table 1). The total number of control versus vinclo- 
zolin lineage differentially expressed genes for all brain re- 
gions combined was 1833 for females and 1693 for males. 
A list of all the genes separated by brain region, sex and 
functional gene categories is presented in Additional file 4: 
Table S2 A-I. The overlap and differences between the 
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Figure 1 A flowchart of the analyses carried out in the study, (a) Differentially expressed (DE) genes are identified by the standard t-test. 
(b) Correlations between the expression profiles of the DE genes are calculated to quantify the level of coexpression. (c) Coexpression network 
analysis is performed to identify coexpressed gene modules. A matrix of correlations between gene expression profiles is first transformed 
through a power function into an adjacency matrix that is further transformed into a topological overlap matrix (TOM). Modules (represented by 
color bars) comprised of highly coexpressed genes are the identified using a dynamic cut-tree algorithm, (d) Association between gene modules 
and phenotypic traits is then accessed, (e) Pathway analysis is then performed on the gene modules of interest to derive regulatory networks for 
mechanism discovery. 



Table 1 Differentially expressed Signature genes and their overlap with modules generated in combined network 
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Signature lists of each brain region for male and female is 
shown in Figure 2. The majority of genes were distinct to 
the different lists in a comparison of the brain regions. The 
one exception was an overlap between the cingulate cortex 
(CngCtx) and olfactory bulb (OlfB) in the female. There- 
fore, each brain region Signature list was distinct from each 
other and between the sexes. 

Analysis of the cellular pathways and processes associ- 
ated with the gene Signature lists for each brain region 
is shown in Figure 2 and Additional file 5. The top 36 
pathways with the greatest combined number of genes 
associated are shown in Additional file 5. An extended 
list of pathways and processes with the associated genes 
from the different modules and tissues is presented in 
Additional file 6: Table S3. Several of the most highly 
represented pathways in the male and female were the 
MAPK signaling pathway, olfactory transduction, neuro- 
active ligand-receptor interactions and axon guidance. The 
Signature list genes distributed relatively evenly across the 
different pathways with no major over-representation iden- 
tified. Most major cellular processes and pathways were 
represented with no major predominance of any individual 
specific regulatory mechanism observed (Additional file 6: 
Table S3). 

A bionetwork cluster analysis was performed on the 
differentially expressed genes in the various brain regions 
as previously described [29,31] to identify gene modules 
and networks with coordinated and interconnected rela- 
tionships (i.e. connectivity) [38] (Figure 1). Initially all the 
differentially expressed genes in the combined brain re- 
gions for male or female were analyzed, termed combined 
networks (Figure 3a). This combined analysis was per- 
formed to potentially identify common gene networks or 



modules similar among all the brain regions that potentially 
correlate with the mate preference behavior parameters. 
The increased number of microarrays and data associated 
with the combined analysis also improves the power of the 
cluster and network analysis. The gene cluster analysis is 
shown and individual modules of genes identified are pre- 
sented in different colors on the axis. The module colors 
represent increasing levels of connectivity [38] with white 
being negligible and red being highest. The combined 
Signature lists provided 4 modules in the female and 9 
gene modules in the male. This can be seen as the blocked 
gene clusters designated as modules of different colors 
(Figure 3a). The number of genes in each module for male 
and female brain regions is shown in Table 1. This com- 
bined network analysis and modules were correlated with 
the mate preference behavior, but no significant correla- 
tions were found using this combined analysis (Additional 
file 7: Table S4). 

Due to the distinct functions of each brain region and 
distinct gene Signature lists, the combined analysis was 
found not to provide the specificity needed to identify the 
behavioral correlations with gene modules. Therefore, a 
more specific network analysis using the individual Signa- 
ture lists for each brain region separately was performed, 
termed separate networks. Each brain region differentially 
expressed gene Signature list was used for separate net- 
work analysis. The cluster analysis identified specific gene 
modules for each brain region from the male and female 
gene Signature lists presented (Figure 3b, c). The modular- 
ity for the specific brain regions was not as strong as the 
combined region analysis. Each separate brain region is 
shown and the gene modules are identified by the differ- 
ent colors. The brain regions had 1-6 different modules 
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Figure 2 Number of overlapped differentially expressed genes with pathways and Signature lists. Number of genes overlapped between 
Signature lists is shown in regular font on grey background; number of affected KEGG pathways overlapped between Signature lists is shown in 
italicized font in white background; only KEGG pathways with 5 or more genes affected are counted. 
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3 Gene Cluster Analysis of Combined Signature Lists Analysis 
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Figure 3 (See legend on next page.) 
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(See figure on previous page.) 

Figure 3 Male and female brain region bionetwork cluster analysis and corresponding gene modules. Topological overlap matrixes of the 
gene co-expression network consisting of genes differentially expressed in F3 generation vinclozolin lineage rat brain as compared to F3 generation 
lineage Control animals. Genes in the rows and columns are sorted by an agglomerative hierarchical clustering algorithm. The different shades of color 
signify the strength of the connections between the nodes (from white signifying not significantly correlated to red signifying highly significantly 
correlated). Modules identified are colored along both column and row and are boxed, (a) Matrixes of combined network for male and female brain 
regions, (b) Matrixes of separate network for female amygdala (F-Amy), cingulate cortex (F-CngCTX), enterorhinal cortex (F-EnCTX), hippocampus 
(F-Hipp), olfactory bulbs (F-OlfB), and preoptic area-anterior hypothalamus (F-POAH). (c) Matrixes of separate network for male amygdala (M-Amy), 
cingulate cortex (M-CngCTX), enterorhinal cortex (M-EnCTX), hippocampus (M-Hipp), olfactory bulbs (M-OlfB), and preoptic area-anterior 
hypothalamus (M-POAH). 



and associated gene networks (Table 1). The same differ- 
entially expressed gene Signature lists were used, but the 
network analysis was from the separate lists (Figure 3b 
and c). All subsequent analyses used the gene modules 
from this region specific network analyses. 

The cluster analysis (Figure 3b and c) for each brain 
region provided modules of genes with coordinated gene 
expression and identified a connectivity index [29-31,33,38] 
for each associated gene. The connectivity index (k.in) for 
each of the differentially expressed genes in each region is 
presented in Additional file 4: Table S2. The top 10% of 
genes with the highest connectivity index for each of the 
gene modules was identified (Additional file 4: Table S2 as 
the genes in bold font). From this combined list of 185 
genes for male and 225 genes for female, a gene sub- 
network analysis was performed. The most highly intercon- 
nected genes in all modules for female and male brain 
regions were used to identify the common direct connec- 
tion interactions between genes in a gene sub-network 
(Figure 4). The female gene sub-network identified angio- 
genesis, growth and apoptosis as predominant cellular pro- 
cesses affected (Figure 4b). The male gene sub-network 
identified apoptosis as a predominant pathway affected 
(Figure 4a). These gene networks identify the common 
connections within the brain regions with the most highly 
interconnected genes differentially expressed between the 
control and vinclozolin F3 generation animals. 

Analysis of the gene networks for each individual brain 
region gene module demonstrated that only males ex- 
hibited direct connection gene sub-networks for gene 
modules (male amygdala and cingulate cortex turquoise 
modules) (Figure 5). This region specific examination of 
gene network modules demonstrated most regions did 
not have direct connection sub-networks, but indirect 
interactions with various pathways and processes. An al- 
ternative analysis used the entire Signature list for each 
brain region to identify region specific gene sub-networks 
(Additional file 8: Figure S2 A-K). How these gene net- 
works may correlate with the alterations in mate prefer- 
ence behavior required a statistical correlation of the gene 
sets with the behaviors (Figure 1). 

In considering the mate preference behaviors, the fe- 
male is the discriminatory sex to choose a mate, while 
the male is non-discriminatory and has phenotypes and 



behaviors to be selected [27,39]. The altered gene expres- 
sion and correlations with behaviors needs to consider 
this in data interpretation. The behavioral parameters 
(Additional file 2: Table SI) for the mate preference analysis 
were statistically correlated to the separate network gene 
modules for the different brain regions, (Additional file 9: 
Table S5). The correlation and the p-values associated with 
the statistical correlation coefficients are presented. All cor- 
relations with a single or multiple principle component 
comparison are presented. Considering a p < 0.05 or correl- 
ation coefficient >0.5 and p = 0.05-0.1 between the gene 
module and behavior demonstrated correlations in four 
female brain regions and six modules with the female 
behavior (Table 2 and Additional file 9: Table S5). A sum- 
mary of the statistically significant correlations and/or 
those with strong correlation coefficients is shown in 
Figure 6. Nearly all the female brain regions had statisti- 
cally significant correlation with the Plexiglas behavioral 
parameter. The female amygdala (F-Amy) had a turquoise 
module with significant correlation with the Walking and 
Still parameters. The turquoise modules of female entorhi- 
nal cortex (F-EnCtx) had a strong correlation with the 
Wire Mesh and Plexiglas behavioral parameters (Figure 6). 

The six male brain regions and associated gene modules 
had a number of statistically significant correlations with 
the mate preference parameters (Table 2 and Additional 
file 9: Table S5). All the male brain regions had statistically 
significant correlation with at least one module and the 
Wire Mesh behavioral parameter. Amy and Hipp also 
had correlations with the behavioral Plexiglas parameter 
(Figure 6). Therefore, at least one gene module in nearly 
all brain regions statistically correlated to the mate pref- 
erence parameters analyzed. These correlations can now 
be considered in regards to the regulatory roles of gene 
networks identified for mate preference behavior alter- 
ations for the female (chooser) versus the male (selected) 
(Figure 5 and Additional file 8: Figure S2). 

The direct connection gene sub-networks for the critical 
male amygdala (M-Amy) and cingulate cortex (M-CngCtx) 
turquoise modules are shown in Figure 5. The Signature list 
for each brain region sub-networks demonstrate distinct 
networks for each region (Additional file 8: Figure S2). 
Since nearly all the brain regions and key modules (Figure 6) 
have a statistically significant correlation with the Wire 



Skinner et al. BMC Genomics 2014, 15:377 
http://www.biomedcentral.com/1471-2164/15/377 



Page 8 of 19 



a Male Direct Connection Gene Sub-Network 




b Female Direct Connection Gene Sub-Network 




Figure 4 Direct connection sub-networks for most highly top 10% connected genes from each module of separate network for male 
(a) or female (b). Only directly connected genes are shown according to their location in the cell (on membrane, in Golgi apparatus, nucleus, 
cytoplasm or outside the cell). Node shapes and color code: oval and circle - protein; diamond - ligand; circle/oval on tripod platform - transcription 
factor; ice cream cone - receptor; crescent - kinase or protein kinase; irregular polygon - phosphatase; red color indicates up-regulated genes, 
blue - down-regulated. Arrows with plus sign show positive regulation/activation, arrows with minus sign - negative regulation/inhibition; grey arrows 
represent regulation, lilac - expression, purple - binding, green - promoter binding, and yellow - protein modification. 
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Figure 5 Direct connection sub-networks for male amygdala turquoise module (a) and male cingulate cortex turquoise module 

Only directly connected genes are shown. Node and arrows shapes and color code is the same as for Figure 4. 



(b). 



Mesh for male or Plexiglas for female mate preference 
behavior parameters, the combined gene sub-network 
(Figure 4) for all male or female brain regions identifies 
a potentially associated molecular control of behavior. 
Alternately, the analysis of separate regions differentially 
expressed gene sets (Signature lists) identified distinct 
gene sub-networks that associate with the different re- 
gions (Additional file 8: Figure S2). These potential gene 
sub-networks correlate and potentially regulate the mate 
choice behavior for the female and selection behavior/ 
phenotype for the male. In addition to the gene networks, 
a correlation of critical cellular pathways in specific brain 
regions and modules (Additional file 5) that are associated 
with the epigenetic transgenerational inheritance of al- 
tered mate preference behavior. 



The epigenetic transgenerational inheritance of the al- 
tered mate preference behavior requires the transmission 
of an altered epigenome in the germline (sperm) [5,13,15]. 
Previously the altered DNA methylation of the F3 gener- 
ation sperm was characterized with 48 differentially DNA 
methylated regions (DMR) being identified in gene pro- 
moters [14]. These vinclozolin induced sperm DMR are 
in part what promotes an altered epigenome in the em- 
bryo and all developing tissues transgenerationally [15]. 
Although any developing tissue (e.g. brain) will have a 
dramatic cascade of epigenetic and genetic steps to 
achieve an adult fully differentiated state [40,41], the 
possibility that some of the original germline epigenetic 
marks (DMR) may persist was investigated. The genes 
associated with the 48 previously identified sperm DMR 
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Table 2 Gene modules highly correlated to different mate preference behavior parameters 



Sex-region Behavior trait 



Wire mesh 



Facial 



Plexiglas 



Still 



Walking 



Module # PC* Correlation p-value Correlation p-value Correlation p-value Correlation p-value Correlation p-value 



F-Amy 


Turquoise 


2 










0.52 


0.028 


0.82 0.012 0.90 


0.001 


F-EnCTX 


Blue 


1 










0.55 


0.066 








Brown 


1 










0.41 


0.037 








Turquoise 


3 


0.81 


0.029 






-0.53 


0.083 






F-OlfB 


Yellow 


2 










0.74 


0.044 






F-POAH 


Turquoise 


1 










0.63 


0.0386 






M-Amy 


Blue 


2 










0.81 


0.008 








Turquoise 


1 


-0.51 


0.022 














M-CngCTX 


Turquoise 


1 


0.54 


0.048 














M-EnCTX 


Turquoise 


1 


-0.58 


0.037 














M-Hipp 


Blue 


1 


0.57 


0.023 






0.77 


0.017 








Turquoise 


1 


0.60 


0.034 














M-OlfB 


Blue 


2 


0.79 


0.012 










0.67 


0.072 




Brown 


1 


0.73 


0.012 
















Green 


1 


0.62 


0.020 
















Red 


1 






0.55 


0.081 












Yellow 


1 


-0.60 


0.033 














M-POAH 


Turquoise 


1 


0.87 


0.002 










0.753 


0.091 



*- number of principal components (PC) used to calculate correlation between modules and behavior. 



were compared with the male and female brain region 
gene Signature lists. The comparison demonstrated the 
majority of the DMR did not correspond to differentially 
expressed genes in various brain regions. Only Rnasel 
in the male Amy, Ig6-2a in the male EnCtx, Parp9 in the 
female CngCtx and Rpl32 in the female OlfB 



overlapped. Interestingly, a copy number variation 
(CNV) in the Famllla site previously identified [14] 
was found in all brain regions identified with the epige- 
nome analysis. This provides a positive control for the 
technology and ability to detect the DMR. Therefore, 
some of the original sperm DMR programmed sites may 
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Figure 6 Gene module correlation with mate preference behavior parameters. Separate network male and female modules highly 
correlated to types of behavior: Wire Mesh (blue bars), Facial Investigation (orange), Plexiglas (green), Walking (red) and Still (black). Bars not 
marked with asterisks have p-value < 0.05; bars marked with one asterisk have correlation coefficient >0.5 and > p-value = 0.1 - 0.05. 
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persist, but the vast majority of brain development and 
epigenetic programming, and potential distal regularity 
role of DMR in epigenetic control regions [25], is dis- 
tinct from the original germline epigenetic marks. 

Discussion 

A systems biology analysis of environmentally induced 
epigenetic transgenerational inheritance of altered mate 
preference behavior was performed to suggest a poten- 
tial role for epigenetics in evolutionary biology. Previous 
research has demonstrated that environmental toxicants 
such as the fungicide vinclozolin can promote a repro- 
gramming of the germline epigenome during fetal gonadal 
sex determination that then transmits altered phenotypes 
and adult onset disease states transgenerationally in the 
absence of future environmental exposure [41]. This is re- 
ferred to as epigenetic transgenerational inheritance [5,15] 
and suggests a role for environmental epigenetics in the 
inheritance of phenotypic variation and disease, independ- 
ent of classic genetic inheritance mechanisms. The basic 
molecular mechanism involved in this non-genetic form 
of inheritance is the ability of environmental factors to in- 
fluence the epigenetic programming of the germline 
[15,19]. The primordial germ cells during migration down 
the genital ridge undergo an erasure of DNA methylation 
that then is initiated to re-methylate at the time of gonadal 
sex determination in a sex-specific manner [19]. An envir- 
onmental toxicant such as vinclozolin appears to alter go- 
nadal development to influence germline DNA 
methylation programming [42] and the differential DNA 
methylation regions (DMR) in the sperm become 
imprinted-like sites that appear to not get erased at 
fertilization so are transmitted to subsequent generations 
and male and female progeny [13-16]. In addition to vin- 
clozolin, a number of other environmental toxicants such 
as the plastic compound bisphenol A (BPA) [16,43], dioxin 
[16,44], methoxycholor [13], phthalates [16], pesticides 
[16], hydrocarbons [16], and DDT [17] have been shown 
to induce transgenerational phenotypes. Other environ- 
mental factors such as nutrition and stress can also pro- 
mote transgenerational phenotypes [45-48]. 

The vinclozolin induced epigenetic transgenerational 
phenotypes previously identified included adult onset rat 
disease after 12 months of age of male infertility, mam- 
mary tumors, prostate disease, kidney disease and im- 
mune abnormalities [20]. Therefore, the mate preference 
analysis was performed prior to adult onset disease to re- 
move the disease as a confounding factor. A brain-behavior 
transgenerational phenotype observed was increased female 
anxiety and decreased male anxiety behaviors [23]. This 
transgenerational anxiety behavior was also examined on a 
molecular level to identify brain region specific changes in 
different gene expression and gene networks associated 
with the behavior [23]. Similar observations were made in 



the analysis of transgenerational stress responses [26]. Inter- 
estingly, in a previous study we found that vinclozolin in- 
duced alterations in mate preference behavior [27] . Females 
from either control or vinclozolin F3 generation lineages 
prefer control lineage males over vinclozolin lineage males, 
whereas no altered mate preference in males was observed 
[27]. While the standard argument would be that the fe- 
males are the discriminating sex and distinguish between 
males on as yet undetermined phenotype characteristic(s), 
it is important to realize that the absence of evidence (in 
the male) does not mean the evidence of absence of male 
involvement since preference is only the first step in a mat- 
ing sequence. That is, under natural circumstances this is 
followed by a mutual decision. Mating in rodents involves 
pheromone and auditory cues produced by both sexes and 
evidence suggests (see below) that it is under such unfet- 
tered conditions that the complementarity of behavior and 
brain are expressed. 

Observations from the current study need to consider 
the effects on the female brain as potentially altering female 
discrimination and preference. The effects on the male 
brain are presumed to be associated with the characteristics 
(e.g. auditory cues and pheromone production) being 
selected. This altered mate preference behavior suggests 
the existence of an environmentally altered epigenetic 
transgenerational inheritance of mate preference behavior 
[15]. The current study was designed to identify the gene 
bionetworks in various male and female brain regions that 
correlate with the behavior of the transgenerational inher- 
itance model. 

A novel gene bionetwork analysis was developed to 
identify gene networks correlated to disease [29]. The 
approach was to use a large number of microarrays to 
identify transcriptomes in specific tissues associated with 
control versus disease individuals in large cohorts. Dif- 
ferentially regulated genes that are coordinately regu- 
lated and having connectivity [38] are clustered in large 
gene sets to identify modules of genes that associate 
with the disease [29,30,33-35] (Figure 1). More recently, 
we have used a similar approach to investigate a normal 
development process to identify gene bionetworks asso- 
ciated with development [31,32]. The primordial follicle 
development in the ovary was investigated to identify a 
network of growth factors and associated signaling systems 
that regulate follicle development [31,32]. This bioinfor- 
matics approach to identify regulatory gene networks was 
used in the current study to correlate brain gene networks 
to mate preference behavior (Figure 1) in an epigenetic 
transgenerational model [27]. The six different brain re- 
gions isolated from F3 generation control and vinclozolin 
lineage females and males were used in a microarray ana- 
lysis to determine the differential gene expression in each 
brain region. The region specific gene sets, "Signature list", 
and associated gene networks were investigated. 
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Analysis of significant pathways and cellular processes 
potentially influenced by the differentially expressed gene 
sets and networks did not identify predominant or over 
represented pathways. One pathway identified that previ- 
ously has been shown to correlate with sexual selection is 
the olfactory transduction pathway [49-52]. However, most 
of the brain regions and specific networks or gene modules 
affected similar pathways with overlap between most. 
Therefore, no specific pathways were identified and most 
major pathways were influenced. A limitation in gene ex- 
pression studies is that individual genes are assigned a spe- 
cific function, but the gene may be involved in functional 
categories. This limitation needs to be considered in any 
gene expression data interpretation, but genome wide tran- 
scriptome analysis for gene sets has been shown to start to 
address this issue [53]. Combined observations suggest that 
the distinct differential expressed gene sets for the different 
brain regions appear to regulate common cellular processes 
and pathways among the brain regions and modules. 

In contrast, analysis of gene networks identified unique 
gene sub-networks and gene modules associated with each 
brain region differentially expressed gene sets. The specific 
inter-connected genes were unique and overall networks of 
connected genes distinct. Therefore, the different functions 
associated with each brain region and associated with the 
altered mate preference behavior were identified. A statis- 
tical correlation of the gene modules for each brain region 



with the different mate preference parameters measured 
identified a number of statistically significant correlations. 
All but one female brain region (Hipp) had statistically 
significant correlations with the Plexiglas behavior param- 
eter. The female behavior directly associates with the dis- 
crimination and mate preference choice. Interestingly, all 
the male brain regions had some modules with statistically 
significant correlations with the Wire Mesh behavior par- 
ameter. The male behavior and/or phenotype (e.g. phero- 
monal production) associates with the selected behaviors 
and characteristics of the non-discriminant sex. Therefore, 
direct correlations with the gene modules, specific brain 
regions and mate preference behavior parameters were 
identified. The specific gene modules and behavioral pa- 
rameters statistically correlated were distinct between the 
sexes and brain regions, but strong correlations of the 
gene networks to the mate preference behavior was estab- 
lished. Interestingly, both the Plexiglas and Wire Mesh are 
indicators of interest and assessment of the stimulus 
animal. 

The gene bionetwork analysis and statistical correlation 
with the mate preference behavior provides insight into 
the molecular basis of how various male and female brain 
regions correlate and in part control the various behav- 
ioral parameters. Observations provide one of the first 
genomic and systems biology analysis of mate preference 
behavior (Figures 1 and 7). The experimental model used 
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Figure 7 Schematic of role of epigenetics in evolution. 
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involved the ability of an environmental compound (vin- 
clozolin) to induce an epigenetic modification of the 
germline (sperm) to promote epigenetic transgenerational 
inheritance of an altered mate preference behavior. This 
altered mate preference behavior was due to a baseline al- 
teration in the epigenomes of all male and female tissues, 
including the brain, which are derived from the epigeneti- 
cally altered germline [25]. The current study used a sys- 
tems biology approach to help elucidate the molecular 
control of this process. 

Darwin [12] considered natural selection and sexual 
selection as distinct processes driving the evolution of 
traits. Natural selection results in traits that are adaptive 
responses to changes in the environment. The resulting 
variation in traits between and within species is shaped 
by differential survivorship. In other words, animals that 
survive are those with traits that are adaptive to their 
environment [54]. Darwin [12] conceived of sexual se- 
lection as arising from aggressive interactions between 
males (male-male competition) and the female's selection 
of a mate (mate choice). Males compete amongst them- 
selves for access to females. Aggression between males can 
have a direct effect on female reproduction by preventing 
other breeding males from having access to females or 
from harming the female. Importantly aggression amongst 
males can also have an indirect effect by inhibiting or sup- 
pressing the normal reproductive physiology of the female 
or even terminate a pregnancy [55]. 

Mate preference, in its simplest form, states that males 
compete for females and females choose between them. 
Although most research has focused on how females 
choose males, male choice of females is also important 
[56-58]. This point cannot be overemphasized. That is, 
in virtually all paradigms published to date, the choosing 
individual is the independent variable and the stimulus 
animal is the dependent variable. Although this study is 
similar to previous studies in that there are restrained 
stimulus animals and freely moving individuals that are 
the investigators, it differs in several important ways. 
First, the 'round robin' testing method insured that all 
males and all females served both as stimulus animals 
and experimental animals; thus, the "preferences" exhib- 
ited reflect both the males and the females. Second, this 
study extends to the molecular level events in the brain 
of the individuals, both of which have exhibited prefer- 
ences, namely brain transgenerational transcriptome al- 
terations that correlated to the opposing sex behaviors. 
In essence we are seeing the outcome of the comple- 
mentary nature of mate choice. This has never been 
shown before. 

Evolution favors reproductive success, and it is in the 
individual's interest to focus on selecting the best mate 
and to avoid mating with the wrong species [59,60]. 
Making the correct choice of a mate has a pronounced 



impact on reproductive success of both partners. Except 
in unusual systems, in nature the mating partners choose 
one another [55,57,60]. Experiments with flies [61], birds 
[62], and rodents [63,64] indicate that individuals who are 
allowed to select, and be selected by, their mate enjoy 
greater reproductive success than force-paired animals. 
This consent is based not only on the internal milieu that 
motivates each individual to seek a partner, but also on 
the satisfactory nature of the phenotypic traits the poten- 
tial mate displays. 

There are a number of sexual selection hypotheses, all 
of which emphasize that females choosing optimally will 
produce young whose viability and survivorship are 
enhanced by the females choice of mate [65]. The most at- 
tractive, and one that takes into account that mating is a 
cooperative act that involves both partners, is the sensory 
exploitation hypothesis [66]. This hypothesis postulates that 
males have evolved calls and/or pheromones to exploit the 
preexisting sensory biases in the female that themselves 
evolved for reasons independent of female choice. Male be- 
havior then has changed to maximize stimulation of the 
female's sensory systems. Therefore, the coordination of the 
complementary signal and receiver, mounting and lordosis, 
coordination of egg and sperm maturation and release is re- 
quired for successful completion for reproduction. These 
complementary processes are evident at all levels of bio- 
logical organization [55,67] and we extend it here to the 
level of the genome and epigenome. 

Conclusions 

The ability of an environmental factor to alter mate pref- 
erence behavior suggests a critical role of environment 
in evolutionary biology. This is distinct from the generally 
accepted role of the environment in natural selection 
where environment is the active factor in the selection of 
an adaptive phenotype, but alternately here it involves the 
induction of phenotypes that can be acted on by natural 
selection. Since the majority of environmental factors can 
not alter DNA sequence or promote mutagenesis [2], an 
additional molecular mechanism to consider involves 
environmental epigenetics [15]. Many environmental 
compounds and factors such as nutrition can modify 
the epigenome to alter phenotypic variation. The role of 
epigenetics in evolutionary biology has been suggested 
previously [5,8-11,15,68], but no significant experimental 
evidence has been provided. The current study demon- 
strates an environmental factor can promote the epigenetic 
transgenerational inheritance of an altered mate preference 
behavior. The epigenetic modification of the germline 
(sperm) has been previously established [13,14,16] and will 
lead to epigenetic alterations in the brain transcriptomes of 
both females and males [23] to alter the mate preference 
behavior [27]. Therefore, the current study provides direct 
experimental evidence for a potential role of environmental 
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epigenetics in evolution by regulating a critical determinant 
such as mate preference on a molecular level (i.e. altered 
gene networks) in specific brain regions in a sex-specific 
manner (Figure 7). Although no direct epigenetic alter- 
ations were examined in the brain, the germline (sperm) 
epigenetic alterations that generate this altered male or fe- 
male brain development have been documented [13,14,16]. 
This molecular mechanism does not suggest genetics will 
not have a critical role in evolutionary biology, but suggests 
environmental epigenetics will be an additional mechanism 
to consider. Epigenetics provides a mechanism for the 
environment to impact phenotypic variation and natural 
selection. Epigenetic and genetic mechanisms will co- 
operate to regulate on a molecular basis evolutionary 
biology. This appears to be a "neo-Lamarckian concept 
to facilitate neo-Darwinian evolution" [40,41]. 

The systems biology approach used in the current study 
links an environmental exposure, epigenetic transgenera- 
tional inheritance and molecular regulation of brain func- 
tion to mate preference and evolutionary biology (Figure 7). 
Epigenetics will have a central role in how environmental 
factors influence how the gene networks emerge to induce 
phenotypic variation. Although genetics is critical for all as- 
pects of biology, epigenetics provides the plasticity to allow 
the environment to alter biological events. This type of 
systems approach to understand complex biological traits, 
such as sexual selection, provides insights into how the 
various components (environment, phenotype and evolu- 
tion) interact in a systems biology manner. 

Methods 

Animal housing protocol 

Male and female rats of the F3 generation of Vinclozolin 
(Vinclozolin-Lineage) and DMSO Control (Control- 
Lineage) Lineages were selected out of litters from un- 
treated F2 generation mothers in Dr. Michael Skinner's 
laboratory at Washington State University according to 
established protocols [13]. Briefly, gestating female FO 
generation Sprague Dawley rats were injected with the 
fungicide vinclozolin (100 mg/kg) daily during fetal go- 
nadal sex determination (E8-E14) and the Fl generation 
were bred to generate the F2 generation and then the F2 
bred to generate the F3 generation [13]. At approximately 
PND 10 (before weaning), each animal was injected with a 
small microchip (AVID Identification system Inc. Norco, 
CA) subcutaneously between the shoulder blades. The 
animals were then shipped to the University of Texas 
from Washington State University on postnatal day (PND) 
22, one day after weaning. Upon arriving at the University 
of Texas, one animal from each Lineage (Control and 
Vinclozolin) was pair-housed (one control and one 
vinclozolin animal) and remained in these dyads through- 
out the duration of the study. Because of the natural vari- 
ation in dates of breeding, there was a 4-day spread of 



birth date of animals in the first cohort but in the second 
cohort, all animals were born on the same day. However, 
all pair-housed animals were no more than one day apart 
in birth age and were paired randomly to prevent an age 
effect on cagemates. 

Each dyad of animals was randomly placed in a six-wide, 
five-high metal housing rack in standard translucent poly- 
carbonate rat cages (46 x 24 x 20.5 cm) with ad libitum 
access to tap water and standard rat chow (Purina rodent 
chow #5LL2 Prolab RMH 1800 diet). The animal room 
was on a 14:10 light/dark schedule. For environmental 
enrichment, a 7 cm diameter PVC pipe was placed in 
each cage. 

Ethics statement 

All experimental protocols for the procedures with rats 
were pre-approved by the Washington State University 
Animal Care and Use Committee (IACUC approval # 
02568-026) and by the University of Texas at Austin 
Animal Care and Use Committee (Public Health Service 
Animal Welfare Assurance Number A4107-01). 

Behavioral testing 

After habituation to the testing arena, each individual 
was tested individually (when used as an experimental 
subject) or in pairs (when used as stimulus animals) with 
all individuals; the order of the testing was rotated dur- 
ing the course of both the male and female trials. All 
tests were conducted during the dark phase of the light 
cycle, beginning at 1200 h, 4 h after the progesterone in- 
jection, in a room illuminated with low levels of red 
light. Before trials, to confirm that females were recep- 
tive, each female was placed with a sexually experienced 
but otherwise experimentally naive male; all females ex- 
hibited robust lordosis (arched back and lifted head pos- 
ture) in response to mounting by the male. 

Partner preference tests consisted of placing an indi- 
vidual (male or female) in the center of a large three- 
chamber glass-testing arena (122 x 46 x 54 cm). At either 
end was a small compartment (28 x 28 x 12.5 cm) con- 
taining the stimulus rats separated by a Wire-mesh barrier 
to allow exchange of olfactory, visual, and tactile cues. The 
area directly in front of the stimulus cage was marked by 
tape. Tests were conducted 2 h after the onset of the dark 
cycle under red-light illumination and lasted 10 min; all 
tests were videotaped for further review and analysis. At 
the end of each test, all animals were removed, and the 
entire testing arena was washed with a household cleaner 
and then wiped down with 70% ethanol to remove scent 
marks and residual odors. All males were tested with both 
types of females as stimulus animals (72 trials), and all 
females were tested with both types of males as stimulus 
animals (72 trials) (Movie SI). 
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The videotaped trials were analyzed by using JWatcher 
vl.O (www.jwatcher.ucla.edu) computer software to quan- 
tify the behavior of each experimental animal. Time spent 
with a stimulus animal was recorded as soon as all four 
paws of the experimental animal crossed over the line of 
tape marking the boundary of that stimulus animal's com- 
partment. As soon as one paw crossed over the tape back 
into the center compartment, the time recorded with the 
experimental animal was stopped. Preference behaviors 
were defined as those directed to the stimulus animal and 
included time spent in contact with the Wire Mesh separ- 
ating the experimental and stimulus animal (Wire Mesh), 
during which the animals often touched noses through 
the Wire Mesh (facial investigation), and contacted the 
Plexiglas surface surrounding the front of the stimulus 
cage; the cumulative total time in these preference be- 
haviors toward each stimulus animal was also calculated 
(Total). Other activity measured included undirected 
walking and sniffing (walking), standing still with minimal 
head movement (still). Videos demonstrating the test can 
be viewed as supporting information (Movie SI). 

Brain processing 

The brain was removed in less than 1 minute and placed 
in crushed ice to chill. The brain was then cut in half in 
the sagittal plane along the midline. In all cases but one 
the right side was blocked and then 6 areas dissected 
(see list below) within 3-5 min. This procedure was done 
on iced tissue. The dissected brain areas were placed in 
chilled Trizol (150 1) in 1.5 ml Eppendorf tubes according 
to manufacturers specifications in each tube. No tissue 
fragment was more than 3 mm but in those instances, but 
where there were multiple fragments the amount of Trizol 
was doubled (approximately). After all animals were 
dissected, the Eppendorf tubes were vortexed (15 sec) 
and then frozen on dry ice. The brain regions collected 
were according to Paxinos & Watson [69]: olfactory 
bulbs (OlfB); cingulate cortex (CngCtx), anterior to POAH 
(Bregma 4.7 to 1.7); preoptic area-anterior hypothal- 
amus (POAH), 4 mm rostral to anterior commissure 
(AC) (Bregma -0.26 to -1.40); amygdaloid nuclei (Amy), 
3 mm caudal to AC (Bregma -2.3 to -3.6); hippocampus 
(Hipp), 6 mm caudal to AC (Bregma -2.12 to -4.52); 
entorhinal cortex (EnCtx) (Bregma -5.60 to -7.80). 

RNA preparation 

Brain area samples from individual rats were homoge- 
nized in 150 ul Trizol and then 600 ul Trizol was added 
to final volume of 750 ml. Samples were stored at -80 
or -20°C until RNA extraction. For microarray analysis, 
from 4 to 6 biological replicas (animals) were prepared 
as above for each brain area Control or Vinclozolin group 
depending on samples availability (Additional file 2: Table 
SIB). A total of 132 (67 Control and 65 Vinclozolin) 



samples/chips were analyzed: (6 brain areas) x (2 Male 
or Female) x (2 Control or Vinclozolin) x (4-6 biological 
replicas). RNA from individual animal brain area was ex- 
tracted from Trizol samples according to standard Trizol 
extraction protocol (Invitrogen, USA) and stored in aque- 
ous solution at -80°C until microarray analysis. 

Microarray analysis 

The microarray analysis was performed by the Genomics 
Core Laboratory, Center for Reproductive Biology, Wash- 
ington State University, Pullman, WA using standard 
Affymetrix reagents and protocol. Briefly, mRNA was 
transcribed into cDNA with random primers, cRNA was 
transcribed, and single-stranded sense DNA was synthe- 
sized which was fragmented and labeled with biotin. Biotin- 
labeled ssDNA was then hybridized to the Rat Gene 1.0 ST 
microarrays containing more than 30,000 transcripts (Affy- 
metrix, Santa Clara, CA, USA). Hybridized chips were 
scanned on Affymetrix Scanner 3000. CEL files containing 
raw data were then pre-processed and analyzed with Partek 
Genomic Suite 6.5 software (Partek Incorporated, St. Louis, 
MO) using an RMA, GC-content adjusted algorithm. Raw 
data pre-processing was performed in 12 groups, one for 
each male or female brain area. Comparison of array sam- 
ple histogram graphs for each group showed if data for all 
chips were similar and appropriate for further analysis 
(Additional file 3: Figure SI). By this criterion, 2 microarray 
samples (not counted in Additional file 4: Table S2B and 
not shown on Additional file 3: Figure SI) were omitted 
from repeated group pre-processing and further analysis. 

The microarray quantitative data involves signals from 
an average of 28 different oligonucleotides (probes) arrayed 
for each transcript and many genes are represented on the 
chip by several transcripts. The hybridization to each probe 
must be consistent to allow a statistically significant quan- 
titative measure of resulting gene expression signal. There- 
fore, the microarray provides an unbiased and highly 
stringent quantitative procedure compared to other proto- 
cols [70]. In contrast, a quantitative PCR procedure uses 
only two oligonucleotides and primer bias is a major factor 
in this type of analysis. Therefore, we did not attempt to 
use PCR based approaches as we feel the microarray ana- 
lysis is more accurate and reproducible without primer bias 
such as PCR based approaches [31]. 

All microarray CEL files from this study have been de- 
posited with the NCBI gene expression and hybridization 
array data repository GEO (GEO series accession number: 
GSE33830) and can be also accessed through www.skinner. 
wsu.edu. For gene annotation, Affymetrix annotation file 
RaGenel_0stvl.na31.rn4.transcript.csv was used. 

Network analysis 

The network analysis was restricted to genes differen- 
tially expressed between the control and the treatment 
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groups based on previously established criteria: (1) fold 
change of group means > 1.2 or < 0.83; (2) Ttest p-value < 
0.05. The union of the differentially expressed genes from 
the different treatments resulted in 1,693 genes for males 
and 1833 for females being identified and used for con- 
structing a weighted gene co-expression network [71,72]. 
Unlike traditional un-weighted gene co-expression net- 
works in which two genes (nodes) are either connected or 
disconnected, the weighted gene co-expression network 
analysis assigns a connection weight to each gene pair 
using soft-thresholding and thus is robust to parameter 
selection. The weighted network analysis begins with a 
matrix of the Pearson correlations between all gene pairs, 
then converts the correlation matrix into an adjacency 
matrix using a power function fix) = xP. The parameter /S 
of the power function is determined in such a way that 
the resulting adjacency matrix (i.e., the weighted co- 
expression network) is approximately scale-free. To meas- 
ure how well a network satisfies a scale-free topology, we 
use the fitting index proposed by Zhang & Horvath [71] 
(i.e., the model fitting index R 2 of the linear model that re- 
gresses log(p(k)) on log(k) where k is connectivity and p(k) 
is the frequency distribution of connectivity). The fitting 
index of a perfect scale-free network is 1. 

To explore the modular structures of the co-expression 
network, the adjacency matrix is further transformed into 
a topological overlap matrix [73]. As the topological over- 
lap between two genes reflects not only their direct inter- 
action, but also their indirect interactions through all the 
other genes in the network. Previous studies [71,73] have 
shown that topological overlap leads to more cohesive and 
biologically meaningful modules. To identify modules of 
highly co-regulated genes, we used average linkage hier- 
archical clustering to group genes based on the topological 
overlap of their connectivity, followed by a dynamic cut- 
tree algorithm to dynamically cut clustering dendrogram 
branches into gene modules [74]. Such networks were 
generated from all combined male or female differentially 
expressed genes (2 combined networks) or from each in- 
dividual male or female brain region Signature lists (12 
separate networks). From one to ten modules were identi- 
fied in combined or separate networks and the module 
size was observed to range from 10 to 780 genes (Table 1). 

To distinguish between modules, each module was 
assigned a unique color identifier, with the remaining, 
poorly connected genes colored grey. The hierarchical 
clustering over the topological overlap matrix (TOM) 
and the identified modules is shown (Figure 1). In this 
type of map, the rows and the columns represent genes 
in a symmetric fashion, and the color intensity repre- 
sents the interaction strength between genes. This TOM 
heatmap highlights that genes in the transcriptional net- 
work fall into distinct network modules, where genes within 
a given module are more interconnected with each other 



(blocks along the diagonal of the matrix) than with genes in 
other modules. Therefore, there are two types of global 
connectivity, adjacency-based one and TO based one. The 
adjacency-based connectivity (kail) is defined as the sum of 
the power-function transformed correlations between the 
gene g and all the other genes in the whole network while 
the TO-based connectivity (to.all) is defined as the sum of 
the topological overlaps between the gene g and all the 
other genes. By default, connectivity used throughout the 
paper refers to TO-based connectivity to.all. 

Gene Co-expression Network Analysis Clarification: 
Gene networks provide a convenient framework for ex- 
ploring the context within which single genes operate. 
Networks are simply graphical models comprised of nodes 
and edges. For gene co-expression networks, an edge 
between two genes may indicate that the corresponding 
expression traits are correlated in a given population of 
interest. Depending on whether the interaction strength 
of two genes is considered, there are two different ap- 
proaches for analyzing gene co-expression networks: 1) 
an unweighted network analysis that involves setting 
hard thresholds on the significance of the interactions, 
and 2) a weighted approach that avoids hard thresholds. 
Weighted gene co-expression networks preserve the 
continuous nature of gene-gene interactions at the tran- 
scriptional level and are robust to parameter selection. 
An important end product from the gene co-expression 
network analysis is a set of gene modules in which mem- 
ber genes are more highly correlated with each other than 
with genes outside a module. Most gene co-expression 
modules are enriched for GO functional annotations and 
are informative for identifying the functional components 
of the network that are associated with disease [75]. 

This gene co-expression network analysis (GCENA) 
has been increasingly used to identify gene sub-networks 
for prioritizing gene targets associated with a variety of 
common human diseases such as cancer and obesity 
[38,76-79]. One important end product of GCENA is 
the construction of gene modules comprised of highly 
interconnected genes. A number of studies have demon- 
strated that co-expression network modules are gener- 
ally enriched for known biological pathways, for genes 
that are linked to common genetic loci and for genes as- 
sociated with disease [33,38,71,75-78,80,81]. In this way, 
one can identify key groups of genes that are perturbed 
by genetic loci that lead to disease, and that define at the 
molecular level disease states. Furthermore, these studies 
have also shown the importance of the hub genes in the 
modules associated with various phenotypes. For example, 
GCENA identified ASPM, a hub gene in the cell cycle 
module, as a molecular target of glioblastoma [78] and 
MGC4504, a hub gene in the unfolded protein response 
module, as a target potentially involved in susceptibility to 
atherosclerosis [77]. 
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J) show gene location in the cell (on membrane, in Golgi apparatus, nucleus, 
cytoplasm or outside the cell). Node shapes and color code: oval and circle - 
protein; diamond - ligand; circle/oval on tripod platform - transcription 
factor; ice cream cone - receptor; crescent - kinase or protein kinase; 
irregular polygon - phosphatase; red color indicates up-regulated genes, 
blue - down-regulated. Arrows with plus sign show positive regulation/ 
activation, arrows with minus sign - negative regulation/inhibition; grey 
arrows represent regulation, lilac - expression, purple - binding, green - 
promoter binding, and yellow - protein modification. 

Additional file 9: Table S5. Correlation between separate network 
modules and behavior trait for F3-Vinlozolin rat brain regions. 



Pathway analysis 

Resulting lists of differentially expressed genes for each 
male or female brain area as well as for each module 
generated in the combined network and some generated 
in separate networks analysis were analyzed for KEGG 
(Kyoto Encyclopedia for Genes and Genome, Kyoto Uni- 
versity, Japan) pathway enrichment using Pathway-Express, 
a web-based tool freely available as part of the Onto-Tools 
(http://vortex.cs.wayne.edu) [82] as well as KEGG website 
'Search Pathway tool (http://www.genome.jp/kegg/tool/ 
search_pathway.html). Global literature analysis of various 
gene lists was performed using Pathway Studio 8.0 soft- 
ware (Ariadne Genomics, Inc., Rockville, MD). 

Additional files 
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